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Abstract 

Rotation-induced splitting of the otherwise degenerate photonic bands is predicted for a two- 
dimensional photonic crystal made of evanescently coupled microcavities. The symmetry-broken 
energy splitting is similar to the Zeeman splitting of atomic levels or electron's (hole's) magnetic 
moment sublevels in an external magnetic field. The orbital motion of photons in periodic photonic 
lattice of microcavities is shown to enhance significantly such Coriolis-Zeeman splitting as compared 
to a solitary microcavity [D.L. Boiko, Optics Express 2, 397 (1998)]. The equation of motion 
suggests that nonstationary rotation induces quantum transitions between photonic states and, 
furthermore, that such transitions will serve as a source of nonstationary gravitational field. 
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The Sagnac effect in a rotating ring cavity, known also as the Coriohs-Zeeman effect for 
photons, emerges as a frequency sphtting of counterpropagating wavesi. Thus, for a ring 
cavity of M wavelengths optical path, the modal shift is MVL. Nowadays, the effect is used 
in commercial He-Ne ring laser gyros of large (M~10^) cavity size. Significant efforts were 
made towards designing miniature- sized solid-state devices based on a high optical gain 
mediunt^. Theoretical investigation has been carried out about the impact of rotation on 
the whispering gallery modes of a microdisc microstructure^. Here, the frequency splitting 
scales with the closed optical path length, while the field polarization is either not relevant 
or assumed to be parallel to the rotation axis f2. 

On the other hand, the Coriolis-Zeeman effect in a cylindrical microwave resonator ro- 
tating along the symmetry axis^ or in an optical Fabry-Perot cavity rotating in the mirror 
plane is, at first sight, independent of the cavity size-. The frequency shift {S+M)VL of 
polarization and transversal modes is set by the spin (±1) and the azimuth mode index 
M. However, the higher is the mode index M, the larger is the size of the cavity that will 
support such mode. By virtue of the complexity of the mode discrimination at high M and 
because of the small frequency splitting of the polarization modes, this effect has not yet 
been verified in experimental measurements. 

In this Letter, I consider the Coriolis-Zeeman effect in coupled microcavities arranged 
in a periodic two-dimensional (2D) lattice in the plane of rotation [Fig{^a)]. On example 
of a square-symmetry lattice, the possibility of enhanced Coriolis-Zeeman splitting, corre- 
sponding to M~1000, is predicted for the low-order photonic modes^. It is caused by the 
photon's orbital motion extended over the large number of lattice cells. The equation of mo- 
tion, which is similar to the Hamiltonian for electrons and holes in magnetic field, suggests 
that nonstationary rotation induces quantum transitions between photonic states and, vice 
verse, that such transitions will generate a nonstationary gravitational field. 

Arrays of evanescently coupled microcavities belong to a particular sub-class of 2D pho- 
tonic crystal (PhC) structures encompassing photonic crystal fibers and arrays of microcav- 
ities. Matrices of vertical cavity surface emitting lasers (VCSELs) are an example of such 
2D photonic crystals''. In such structures, only a small transversal component k_i_ of the 
propagation vector k undergoes Bragg reflections in the plane of periodic lattice [FigiU^a)]. 
The structures employ lattices of periods significantly exceeding the optical wavelength. 
They are typically realized by the mirror reflectivity patterning in a broad-area microcavity. 
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As a model system for such PhC structures, a Fabry-Perot cavity with patterned mirror 
reflectivity is considered here. The cavity length is one wavelength, l-^=\jn with n being 
the refractive index in the cavity. The reflectivity Ri{x,y) of the one cavity mirror {e.g., of 
the top mirror) is modulated in two directions parallel to its plan. The reflectivity pattern 
Ri= expi(f{x,y) consists of pixels with the relative phase shift !f{x,y)=Aip separated by a 
grid of !f{x,y)=0. The mirror is thus perfectly reflecting (|-Ri| =1) and the phase modula- 
tion pattern (p{x, y) defines the structure of the cavity modes. The period of the reflectivity 
pattern A is of a few micron pitch (A^Z^). The pattern is characterized by a fill factor 
FF that is the ratio between the area of the pixel and that of the unit cell. Like in typical 
VCSEL arrays, the phase contrast of reflectivity pattern is small (|Ay9| ^10"^). The second 
mirror of uniform reflectivity (i?2=l) has no impact on the cavity modes. 

The analysis is carried out here using an equivalent, unfolded cavity representation^. 
Multiple reflections at the cavity mirrors effectively translate the cavity into a structure that 
is periodic along the cavity axis (the 2;-axis). The unfolded PhC is thus three-dimensional and 
it can be analyzed in terms of a modal expansion on orthogonal plane waves (OPWs). The 
Coriolis-Zeeman effect in such photonic crystal is considered here using a frame of reference, 
which rotates together with the crystal. Such noninertial rest frame is characterized by a 
metric tensor gik, the off-diagonal space-time components goa of which are dependent on 
the angular rotation speed OS. However, in the unfolded-cavity representation, go. differs 
from the diagonal Minkowski tensor, even in the absence of rotation. Thus goay^O at the 
subsequent mirror reflections since the reflection operator a=IC2 includes rotation by vr 
(about the z axis) followed by coordinate inversion [FiglT](a)]. Due to the patterned mirror 
reflectivity, the equivalent unfolded PhC is of periodically varying " noninertiality" in the xy 
plane^. In the approximation of the first-order terms Vlr/c and Ay9, the metric tensor gik 
has the following nonzero components: 



where goa=gao (« = 1, 2, 3), the space-time coordinates are indexed according to the intervals 
dx^=cdt, dx^=dx, dx'^=dy and dx^=dz; twice repeated Greek indexes indicate summation. 
The first term in goa accounts for rotation of the coordinate systemi. The second term 
accounts for the multiple cavity roundtrips along the 2;-axis and reflections at the cavity 



goo 



-fl'll = -5'22 = -fi'33 = l 



(1) 
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mirrors. The z-period of the unfolded crystal is thus 2/^. The metric tensor ([T]) is validated 
by inspecting the system Hamiltonian [Eq.([7])] for the case of ip=0 (rotating FP cavity^) or 
Q=0 (PhC in an inertia framed). 

In the approximation ([T]), the coordinate space is Euclidean, with the metric tensor 
7q/3= — ga/3 + goaQoisI Qoo being the Kronecker delta dap- Proceeding in a standard manner-, 
the covariant Maxwell's equations with metric ([1]) are converted to the usual form in terms 
of noncovariant field vectors B, H, D and E that assume the constitutive equations 



D=eE + Hxg, B=/iH + gxE, 



fi X r . c , , 



C UJ 



(2) 
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where z is the unit vector along z-axis direction and the components of the vector g are 

ga= — goa/goo- 

Maxwell' equations in photonic crystal ([2]) are solved here by separating fast oscillations 
in the z-axis direction and slow lateral field oscillations in the xy plane: 



Ea 
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(3) 



where n=y/e]I and Z=^y fi/e are the refractive index and impedance in the cavity. The 
gauge transformation is introduced here through the operator 



£al3—Sa/3 ( 1 
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4A;2 dxydx^ 
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2kl dxadxfs 
Xfy d 



(4) 



where the terms ^k\/k1 are taken into account. Such separation of variables is valid in con- 



ditions of the paraxial approximation 



dxf, 



■Cl) and of the low contrast 



of reflectivity pattern {\Aip\ <^1). The two-component vector ■j/'(r_|_)= ^ j in the x-y plane 
is a slowly-varying function of coordinates. It defines the spatial patterns of the six electro- 
magnetic field components ([3]) and it is considered here as the photonic state wave function. 
Its squared modulus \il^{rj_)f yields the intensity pattern of the main polarization compo- 
nent in ([3]). For Q=0, Eqs.([3])-(jl]) are in agreement with the results obtained for the Gaussian 
beam^. 
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In photonic crystals, the wave function is a Bloch wave propagating in the xy 

plane^, 



1^1 



(5) 
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where ^ /^^^^ u*,j^^Ugkj^'^^r±='5g'^- The longitudinal part in ([3]) [the term e 
also a Bloch function. Within the ^-period of the lattice, it has a small modulation depth 
(|^(jk|)2i =2r /-I l^ekl dz^A(p, which is set by an effective phase shift a^k at each reflection 
of the patterned mirror. The general form of such periodic function rjq^ is 



l+?7gk=exp < iaqk^ 



2L 



(6) 



where 9{z)= S{()d( is the unit step function. Note that ?7gk(-2^) is the odd function and 
{dr)q]i/dz)^i —0 by virtue of the small contrast of the reflectivity pattern. 

By operating with cs/jq^"^^ and [from (jlj)] on Maxwell' equations for the curl of E and 
H, substituting the gauge (I3l) and averaging over the 2;-period of the lattice, the Maxwell' 
equations are converted into the same form of a Hamiltonian eigenproblem with respect to 
the photonic state wave function 'j/'(/k(r±) 

^ (7) 



rrioc 



p1 

2mn 
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- ( r± X P± + ^5*2 ) V'gk = h^gklpqk, 



° q) in the 



where mQ=nhkz/c and Sz=izx is the spin operator that reads (5*^) 0/3=^603^ 
basis of the two-component vector functions 1/'= (^^l^- Solutions of Eq.([7]) define the slowly- 
varying components of photonic modes ([3]). The difference between the exact equation for 
(l+'7gk)''/'(jk and its 2;-period average [Eq.(I71)] yields the equation for the periodic part of the 
fast longitudinal component: 



dz 



iaqi, (1 + r^gk) 



1 

2L 



(8) 



where agk= ("^gkl V I'^gk)- The solution of ([8]) is given by ([6]) provided that the eigen functions 
i/'gk of the Hamiltonian ([7]) are known. 

In (l7l), the first term is due to the paraxial propagation along the z axis. The second and 
third terms are the in-plane kinetic energy and the periodic crystal potential, respectively. 
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The last term in ([7]) is a perturbation induced by the Coriohs force. For (f=0 and Q=0, 
within the accuracy of the time variable, Eq.([7]) is just the paraxial wave equation. For ip=0, 
Eq.© yields the effective refractive index ^(k^+j^) of a circularly-polarized paraxial wave 

nes = n + gT± , (9) 

un 

where t=| defines the propagation direction, +/— sign is for the left /right handed po- 
larization. Eq.Q agrees with previously reported expressions for the axial nonreciprocity^ 
(second term) and circular birefringence^ (third term) induced by the Coriolis force for 
photons. Finally, the unperturbed (i7=0) Hamiltonian has been verified by experimental 
measurements in VCSEL array PhC heterostructuresii. These justify the approximation 
of the Hamiltonian. 

Analytical similarities between the effective single-electron Hamiltonian in a semiconduc- 
tor subjected to an external magnetic field and Eq.(I7j) allow the correspondence between 
the periodic crystal potential V and phase pattern ip iy^—-^v)i and the vector poten- 
tials A=iHxr and g_L=irixr (^A^^^gx). Photons in photonic crystal subjected to a 
nonpermanent gravitatinoal field exhibit thus a similar behaviour with electrons (holes) in 
a magnetic field. Accordingly, the impact of rotation on the envelope function and peri- 
odic part of the photonic Bloch wave ([5]) is different. Like in the case of electrons,-^ the 
components of velocity operator \^ = ^[y^,H] do not commute (['Oi^,!)^] = — :^^^VLz^ where 
v±=~-~ "^2"^ )- However, the second-order terms in the Hamiltonian ([7]) [and, re- 

spectively, in ([1])] are needed to define whether the Landau-like quantization is possible for 
photonic envelope wave functions. In the rest of the Letter, the impact of rotation on the 
periodic part of Bloch functions is examined in details on example of a square-lattice PhC. 

FiglT](b) shows the typical band structure of a square-lattice PhC, which is calculated 
along the high symmetry lines A-Z-T in the Brillouin zone (BZ), using the OPW expansion 
in unperturbed (i7=0) Hamiltonian. By virtue of the square lattice symmetry, all states 
are degenerate by the photon's spin {e.g., the doubly degenerate states T5 or Tg). Angular 
rotation removes the degeneracy of such states and splits their energies on 2^ [Figl2](a), 
top panel]. However, there are states, like the degenerate states Ti, T2, T3 and T4, of the 
four-fold degeneracy, which is caused by the orbital symmetry of the Bloch functions. For 
such states, there is an important angular momentum contribution -^hVtzLz to the energy 
shift. 
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The Coriolis-Zeeman splitting of these states is analyzed here using the first-order kp 
expansion in the T point [k = ^, k^)] of the BZ. The expansion basis is deduced from 
the empty lattice test. Four scalar plane waves exp(±^±^), which form the first photonic 
band of empty lattice, originate from the nearest equivalent T points of reciprocal lattice 
and provide representation that is reducible under the C^^ point group (the symmetry group 
of k). Their symmetrized combinations of Ti, T4 and T5 representations are indicated in 
the second column of Table [T] with the capitals letters corresponding to the main term of 
Taylor expansion on parameter {e.g., jS") =^ cos ^ cos =^ sin ^ cos ^). The 

photon's spin transforms as the two-dimensional representation T5. Therefore, a reduction 
of the direct product Tj ® T5 results in the eight symmetry adapted photonic harmonics of 
T1-T5 and Tg representations that constitute a suitable kp-expansion basis in the low-order 
photonic bands. These states are well separated energetically from the other states in the 
T point [Fig. 1(b)]. 

For a general state \il^qk) =6*^^^"-^ q |Tj) at k_i_ measured from the T point of the BZ, 
the 8x8 kp Hamiltonian for the coefficients q is of the block- diagonal form 



H 
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where 
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in the basis of functions ^in.^^ny) , ^(TiTiTa) , ±|(T3±iT4) , ^iCTsxiiTsy)^. The up- 
per (lower) sign refers to the top (bottom) 4x4 block. In (fTOj) . k±=kx ± iky, P= 
^(5*1— z^^|zX)=-^^ is the interband matrix element of p that defines the band mix- 
ing. The eigen solutions of (ITU]) are plotted in Fig. 1(b). Good agreement with the band 
structure calculated from Eq.([7]) justifies the first-order kp approximation flTUl) of the system 
Hamiltonian. 

The matrix Hq in (fTO!) is a perturbation induced by rotation. It accounts for the Coriolis- 
Zeeman energy shift, which is of the opposite sign for the left (upper 4x4 block) and right 
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(lower block) handed polarization states of the same orbital symmetry. The rotation lifts 
degeneracy between such spin states, producing the frequency splitting Aus = 2^, as indi- 
cated in FigEt^a) [top panel] on example of the T5 (T5) states. In ( fTOl) . the parameter M= 
M+H-M. accounts for the orbital part of wave functions. The orbital contribution Aujl='2^^ 
to the Coriolis-Zeeman energy shift ±^hAu!L±^hAus in the T1-T4 states is shown in the 
bottom panel of Figl2](a) {Aus is not visible at the scale of Aujl). The energy shift ±lhAujL 
[the term — ^rxp in Eq.©] is evaluated here using the relationship ^mn=^ ^^T^ ■ Note 
that the procedure to evaluate the matrix elements of hLz=rj_x-pj_ in free spaco^ and in 
periodic latticea^ is different. Here, the f-sum rule ^— =^-| — ^ V ^"^g"* — implies that 

f ' J niap mo mg En-Em ^ 

M± = tUi^-^) and 



Aul = —r 



rriQ rriQ 



A.. = 2§. (11) 



For an array of square pixels defining the microcavities, 

1 

It follows that reducing the effective mass 



hmocFFAip 



sin ttV FF / 1 _i_ sin vr v FF ^ 



/FF ^ nVFF 



(via the lattice pitch A, fill factor FF and contrast A(p), one can enhance the Coriolis- 
Zeeman splitting and achieve M of more than 10^. [FigEKb) shows the ratio ^ = left 
axis]. The enhancement is caused by the weak localization of photonic wave functions to 

the lattice sites. The intraband matrix element (gk|r^|gk)= — — in the T1-T4 bands 

indicates that the wave functions (and, hence, the photon's orbital motion) spread over a 
large PhC crystal domain [Figl2](b) right axis], such that ^ = '"^ff^ [l+sinc^7rVFF]-^ 
In the bands of T5 symmetry, the intraband matrix elements of are nonzero as well, 
however, the orbital contribution to the Coriolis-Zeeman energy vanishes, in accordance 
with the group theory selection rules. With present experimental techniques, the predicted 
frequency splitting [Figl2](a)] can be validated by direct measurements. 

The analogy between Eq.(l7l) and electron's (hole's) magnetic moment Hamiltonian sug- 
gests that a nonstationary field g(t) will stimulate transitions between the photonic states, 
and that such transitions will serve as a source for nonpermanent gravitational field g(t). 

The author is grateful to Eli Kapon, Marc-Andre Dupertuis and Edoardo Charbon, as 
well as to the Fonds National Suisse de la Recherche Scientifique. 
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Figures 




FIG. 1: (a): Schematic of the array of coupled micro cavities, (b): Band structure of square- 
lattice PhC calculated from ([7]) [black curves] and using 8x8 k-p approximation (jlOp [gray curves]. 
The inset shows the BZ and location of the high-symmetry triangle AZT. The parameters are 
A=960nm, n=3.53, A=4/um, FF=0.65 and Av9=0.02. 
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FIG. 2: (a): Coriolis-Zeeman splitting of the T5 (T'5) (top panel) and T1-T4 (bottom panel) 
bands as a function of rotation rate Q. A = 4^m and Aip = 10~^. (b): Relative splitting ^'^q'^ 
(left axis) and matrix element (Ti|r^|Ti)2 (right axis) as a function of the lattice contrast Aip. 
The lattice constant A is 4 (black curves) and 6^m (gray curves), n = 3.53. Other parameters are 
given in the caption of FigHJ 
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Tables 



TABLE I: Basis functions (scalars) and photonic harmonics (vectors) of irreducible representations 
of the group of k at the T point of the BZ (C4^ point group) 



T, T, X T, 

Ti S T5 S±,Sy 

T2 XY{X^-Y^) T'l XY{X^ -Y'^)±,XY{X^ ~Y^)y 

T3 X2-y2 {X^ -Y^)i,{X^ -Y^)y 

Ti XY T'5 XYyi,XYy 

n iX,iY T1 + T2 + T3 + T4: 

Ti ^{XSc + Yy) 

T2 -^(Y5c-Xy) 



T3 ^ (Xx - Yy) 

T4 -j={Y± + Xy) 



12 



